Stromal Signals Dominate Gene Expression Signature Scores That Aim to Describe Cancer Cell-intrinsic Stemness or Mesenchymality Characteristics.

Epithelial-to-mesenchymal transition (EMT) in cancer cells confers migratory abilities, a crucial aspect in the metastasis of tumors that frequently leads to death. In multiple studies, authors proposed gene expression signatures for EMT, stemness, or mesenchymality of tumors based on bulk tumor expression profiling. However, recent studies suggested that noncancerous cells from the microenvironment or macroenvironment heavily influence such signature profiles. Here, we strengthen these findings by investigating 11 published and frequently referenced gene expression signatures that were proposed to describe EMT-related (EMT, mesenchymal, or stemness) characteristics in various cancer types. By analyses of bulk, single-cell, and pseudobulk expression data, we show that the cell type composition of a tumor sample frequently dominates scores of these EMT-related signatures. A comprehensive, integrated analysis of bulk RNA sequencing (RNA-seq) and single-cell RNA-seq data shows that stromal cells, most often fibroblasts, are the main drivers of EMT-related signature scores. We call attention to the risk of false conclusions about tumor properties when interpreting EMT-related signatures, especially in a clinical setting: high patient scores of EMT-related signatures or calls of "stemness subtypes" often result from low cancer cell content in tumor biopsies rather than cancer cell-specific stemness or mesenchymal/EMT characteristics. SIGNIFICANCE: Cancer self-renewal and migratory abilities are often characterized via gene module expression profiles, also called EMT or stemness gene expression signatures. Using published clinical tumor samples, cancer cell lines, and single cancer cells, we highlight the dominating influence of noncancer cells in low cancer cell content biopsies on their scores. We caution on their application for low cancer cell content clinical cancer samples with the intent to assign such characteristics or subtypes.

Immune cell type signature discovery and random forest classification for analysis of single cell gene expression datasets.

Background: Robust immune cell gene expression signatures are central to the analysis of single cell studies. Nearly all known sets of immune cell signatures have been derived by making use of only single gene expression datasets. Utilizing the power of multiple integrated datasets could lead to high-quality immune cell signatures which could be used as superior inputs to machine learning-based cell type classification approaches. Results: We established a novel workflow for the discovery of immune cell type signatures based primarily on gene-versus-gene expression similarity. It leverages multiple datasets, here seven single cell expression datasets from six different cancer types and resulted in eleven immune cell type-specific gene expression signatures. We used these to train random forest classifiers for immune cell type assignment for single-cell RNA-seq datasets. We obtained similar or better prediction results compared to commonly used methods for cell type assignment in independent benchmarking datasets. Our gene signature set yields higher prediction scores than other published immune cell type gene sets in random forest-based cell type classification. We further demonstrate how our approach helps to avoid bias in downstream statistical analyses by re-analysis of a published IFN stimulation experiment. Discussion and conclusion: We demonstrated the quality of our immune cell signatures and their strong performance in a random forest-based cell typing approach. We argue that classifying cells based on our comparably slim sets of genes accompanied by a random forest-based approach not only matches or outperforms widely used published approaches. It also facilitates unbiased downstream statistical analyses of differential gene expression between cell types for significantly more genes compared to previous cell classification algorithms.

Gene expression-based prediction of pazopanib efficacy in sarcoma.

BACKGROUND: The multi-receptor tyrosine kinase inhibitor pazopanib is approved for the treatment of advanced soft-tissue sarcoma and has also shown activity in other sarcoma subtypes. However, its clinical efficacy is highly variable, and no reliable predictors exist to select patients who are likely to benefit from this drug. PATIENTS AND METHODS: We analysed the molecular profiles and clinical outcomes of patients with pazopanib-treated sarcoma enrolled in a prospective observational study by the German Cancer Consortium, DKTK MASTER, that employs whole-genome/exome sequencing and transcriptome sequencing to inform the care of young adults with advanced cancer across histology and patients with rare cancers. RESULTS: Among 109 patients with available whole-genome/exome sequencing data, there was no correlation between clinical parameters, specific genetic alterations or mutational signatures and clinical outcome. In contrast, the analysis of a subcohort of 62 patients who underwent molecular analysis before pazopanib treatment and had transcriptome sequencing data available showed that mRNA levels of NTRK3 (hazard ratio [HR] = 0.53, p = 0.021), IGF1R (HR = 1.82, p = 0.027) and KDR (HR = 0.50, p = 0.011) were independently associated with progression-free survival (PFS). Based on the expression of these receptor tyrosine kinase genes, i.e. the features NTRK3-high, IGF1R-low and KDR-high, we developed a pazopanib efficacy predictor that stratified patients into three groups with significantly different PFS (p < 0.0001). Application of the pazopanib efficacy predictor to an independent cohort of patients with pazopanib-treated sarcoma from DKTK MASTER (n = 43) confirmed its potential to separate patient groups with significantly different PFS (p = 0.02), whereas no such association was observed in patients with sarcoma from DKTK MASTER (n = 97) or The Cancer Genome Atlas sarcoma cohort (n = 256) who were not treated with pazopanib. CONCLUSION: A score based on the combined expression of NTRK3, IGF1R and KDR allows the identification of patients with sarcoma and with good, intermediate and poor outcome following pazopanib therapy and warrants prospective investigation as a predictive tool to optimise the use of this drug in the clinic.

CONET: copy number event tree model of evolutionary tumor history for single-cell data.

Copy number alterations constitute important phenomena in tumor evolution. Whole genome single-cell sequencing gives insight into copy number profiles of individual cells, but is highly noisy. Here, we propose CONET, a probabilistic model for joint inference of the evolutionary tree on copy number events and copy number calling. CONET employs an efficient, regularized MCMC procedure to search the space of possible model structures and parameters. We introduce a range of model priors and penalties for efficient regularization. CONET reveals copy number evolution in two breast cancer samples, and outperforms other methods in tree reconstruction, breakpoint identification and copy number calling.

Deconvolution of sarcoma methylomes reveals varying degrees of immune cell infiltrates with association to genomic aberrations.

BACKGROUND: Soft-tissue sarcomas (STS) are a heterogeneous group of mesenchymal tumors for which response to immunotherapies is not well established. Therefore, it is important to risk-stratify and identify STS patients who will most likely benefit from these treatments. RESULTS: To reveal shared and distinct methylation signatures present in STS, we performed unsupervised deconvolution of DNA methylation data from the TCGA sarcoma and an independent validation cohort. We showed that leiomyosarcoma can be subclassified into three distinct methylation groups. More importantly, we identified a component associated with tumor-infiltrating leukocytes, which suggests varying degrees of immune cell infiltration in STS subtypes and an association with prognosis. We further investigated the genomic alterations that may influence tumor infiltration by leukocytes including RB1 loss in undifferentiated pleomorphic sarcomas and ELK3 amplification in dedifferentiated liposarcomas. CONCLUSIONS: In summary, we have leveraged unsupervised methylation-based deconvolution to characterize the immune compartment and molecularly stratify subtypes in STS, which may benefit precision medicine in the future.